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YUXI HU AND YACHUN LI 

Abstract. We study a Neumann type initial-boundary value problem for strongly 
degenerate parabolic-hyperbolic equations under the nonlinearity-diffusivity con- 
dition. We suggest a notion of entropy solution for this problem and prove its 
uniqueness. The existence of entropy solutions is also discussed. 



1. Introduction 

In this paper, we consider the following initial-boundary value problem: 

Ut + V -iiu) = AB{u), (x,t) e X (0,T), 
u\t=o = Mo(x), X G (1.1) 
(f (n) - VB{u)) ■ n(x, t) = 0, (x, t) e dQ x (0, T), 

where fl is an open and bounded rectangle of , T > and n(x, t) denotes the 
outward unit normal vector to dQ at (x, t). tio(x) : — )■ M is a bounded function 
and B{u) is a non-decreasing function in C^(M). The flux function i{u) is assumed 
to be regular from R to M^. 

Partial differential equations of the above form arise in mathematical models of 
many physical situation, such as two phase flows in porous media, sedimentation- 
consolidation process, etc (see [1], [H], [T5]). 

Well-posedness of general Dirichlet boundary value problem as well as Cauchy 
problem for both isotropic and anisotropic degenerate parabolic-hyperbolic equa- 
tions have been well studied by [12], [16], [19], [21], [22]. However, the Dirichlet 
boundary condition may not always provide the most natural setting for various 
physical problems, see [6J, For example, for equations modeling sedimentation- 
consolidation processes, it's usually appropriate to use kinematic "flux- type" or 
"wall" boundary conditions rather than the Dirichlet boundary condition (see [S]). 
In one-dimensional case, specific models (sedimentation-consolidation model) have 
been considered by Biirger, Evje, Frid and Karlsen [6] [7]. In [7j, well-posedness 
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has been established under the "flux type" boundary condition using BV approach. 
It turned out that these boundary conditions are satisfied on the boundary in an 
a.e. pointwise sense, which is difficult to be extended to the multidimensional case 
due to lack of regularity on diffusion term B{u). On the other hand, although the 
BV approach plays an important role in the study of degenerate parabolic equa- 
tions, it does have many restrictions, see [B]. Based on such considerations. Burger, 
Frid and Karlsen [6] considered a free boundary problem under the framework of 
divergence-measure fields which was first considerd by Anzellotti [2] and reformu- 
lated by Chen and Frid [10] [11]. It should be noted that the boundary conditions 
proposed in [6] is too weak to ensure the uniqueness of entropy solutions. It moti- 
vated us to impose some stronger boundary condition but still under the framework 
of divergence-measure fields. Moreover, to our knowledge, there are few results refer- 
ring to Neumann boundary conditions for strongly degenerate parabolic equations 
in multidimensional case. We point out that for a kind of special weakly degenerate 
parabolic-hyperbolic equation with B'{u) > for u ^ and B'{0) = , Anderson 
[3] has proved local well-posedness with mixed-boundary value condition. However, 
since we can not expect any regularity results for strongly degenerate parabolic- 
hyperbolic equation, the method in |3] can not be used here. 

When B[u) is strictly increasing, the above equation is of parabolic type. On this 
occasion, existence and uniqueness of weak solutions are well-known, see [THj. In the 
case where B' = 0, problem fll.ip is reduced into a nonlinear hyperbolic conservation 
laws with a zero-fiux boundary condition. This kind of problem has been recently 
studied by Biirger, Frid and Karlsen where they introduced a new formulation 
of entropy solutions and proved its well-posedness. A main feature in this new 
formulation lies in the existence of strong traces of any L°° entropy solutions, which 
is based on the result of Vasseur [2l]. It should be noted that there is a mistake in 
the proof of uniqueness in [5j which was solved by Hu and Li [13] by using a different 
approach. 

This paper is mainly devoted to giving a proper formulation of entropy solutions 
for Neumann boundary value problem f 1 1.11) and proving its uniqueness only on rec- 
tangle domain. The main difficulty is to explain the meaning of boundary condition 
(f(7i) — VB{u)) ■ n(x, t) = 0. Since Uo{x) is only supposed to be bounded in do- 
main Q, one can not expect the solution m(x, t) belonging to a more regular space 
than L°° space. Therefore, the meaning of f(tt) — VB{u) on the boundary is not 
clear. Fortunately, by theory of divergence-measure field introduced by Chen and 
Frid [lOj, ^IJ, the vector field (m, f (u) - VB{u)) e DM^{{0,T) x n) has normal 
trace on the boundary. Thus, the boundary condition fll.ll) o can be understood as 
that the normal trace of f{u) — VB{u) is zero on the boundary. However, as far 
as we know, there are no results concerning uniqueness of entropy solutions under 
such boundary condition. Therefore, we impose some stronger conditions on the 
boundary which in turn need some stronger regularity for solutions, that is, the 
existence of strong traces on the boundary. For this purpose, we impose so-called 
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nonlinearity-diffusivity conditions on the flux function i{u) and diffusive function 
B{u). Tliis condition is used to ensure tlie existence of boundary traces via tfie 
result of Kwon [17] as well as the compactness via the result of Lions, Perthame and 
Tadmor ^U\. Our boundary condition is reasonable since it can be derived though 
an entropy inequality, and thus can be understood as an entropy boundary condi- 
tion. We also point out that if the solution u is regular, our boundary condition is 
consistent with classical boundary conditions. 

The reason why we consider rectangle domains is as follows: The boundary condi- 
tion we imposed is still too weak to ensure uniqueness of entropy solutions for general 
domains. The special structure of rectangle domains allows us to understand the 
boundary condition in an a.e. pointwise sense in some part of the boundary dQ. 
Therefore, the boundary terms can be dealt with and uniqueness can be derived. 

The existence of entropy solutions is also discussed. For some special flux function 
f and diffusion function B, we can show that the entropy solution exists by using the 
method of vanishing viscosity. More specifically, we assume that there is a critical 
value Uc such that B'{^) = for ^ < Uc and B'{^) > for ^ > Uc and vanishes 
on ^ > Uc- Here we do not claim that the flux function and the diffusion function 
satisfying these conditions are appropriate for any real physical models. However, 
they do satisfy our assumptions for well-posedness of problem (11. ip . We also note 
that such chosen B{^) is natural for sedimentation-consolidation processes, see [6], 

We organize our paper as follows. In section 2, we state some technical assump- 
tions on the flux f (w) and introduce the concept of domain with Lipschitz deformable 
boundary as well as the notations of boundary-layer sequence and normal traces of 
divergence-measure field. Meanwhile, we present a strong trace results mainly from 
Kwon [17] . In section 3, we give the definition of entropy solutions to problem (11. ip 
and prove its uniqueness. In section 4, we discuss the existence of entropy solutions. 

2. Preliminaries 

We need some assumptions on the flux f('u) as in [5]. 

Assumption 2.1. Assume that uq G [0, M] for fixed M > and i{u) depends 
smoothly on u for u G [0, M]. Moreover, we require that 

f(0)=0, f(M) = 0. (2.1) 

Assumption 2.2. Assume that the flux function i{u) and the diffusion function 
B{u) satisfy the nonlinearity-diffusivity condition, that is, V(r, ^) G M x ]R^,r^ + 

£({m G [0,M] |t + ^ ■£'(«) = and B\u)\i\^ = Q}) = Q , (2.2) 

where L denotes the one- dimensional Lebesgue measure. 

Remark 2.1. As in [5], Assumption \2.1\ is used to ensure L°° hounds on the solu- 
tions. 
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Remark 2.2. The nonlinearity-diffusivity condition implies that there are no in- 
terval of u on which both f{u) is affine and B{u) is degenerate. This condition is 
mainly developed from [131 120] and is always used to the theory of compensated com- 
pactness [23]. Moveover, as is shown in Lemma \2.1\ the condition fl2.2p implies the 
existence of strong traces of entropy solutions. 

2.1. Deformable Lipschitz boundary and Strong traces. Now we introduce 
the concept of deformable Lipschitz boundary [lOj and the notion of strong traces. 

Definition 2.1. We say that dQ is a deformable Lipschitz boundary provided that 
the following conditions hold: 

(1) For all X G dVt there exists a number r > and a Lipschitz map h : M^^^ — )• R 
such that, after rotating and relabeling coordinates if necessary, 

n n Q(x, r) = {y G : h{y,, ■ ■ ■ , vn-i) < Vn} H Q(x, r), 

where Q(x, r) = {y G : \xi — yi\ < r, i = 1, ■ ■ ■ , N}. We denote by h the map 

ivi,--- ,yN-i)=-y^{y,h{y))- _ 

(2) There exists a map ^ : dVL x [0, 1] Vt such that is a homeomorphism, that 
is, bi-Lipschitz over its image with ^'(a;, 0) = u for all u G dQ. 

The map \1/ is called a Lipschitz deformation of the boundary dQ. We denote 
^s(a;) = ^(a;,s) and dQg = ^s(5fi). We also denote by Qg the bounded open 
set whose boundary is dQs- 

Moreover, the Lipschitz deformation \1/ is said to be regular if 

hmV^, o/i = V/i inLl^iB), (2.3) 

s— >-0+ 

where B denotes the greatest open set such that h{B) C. dfl. 
The meaning of strong trace is stated in the following sense: 

Definition 2.2. Let M. C M^+^ have a regular deformable Lipschitz boundary. 
We say that a given function u G L°°{M.) possesses a strong trace on dAi 
if G L°°{dM.) have the property that for every regular (with respect to dAi) 
Lipschitz deformation \E' and every compact set K CC dM., 

essWm I |m(^(s,x)) -ii"(x)|d?/^(x) = 0, (2.4) 
Jk 

where l-L^ is the N -dimensional Hausdorff measure. 

From now on in this paper, it is always understood that VL C M.^ is an open 
bounded rectangle and thus has deformable Lipschitz boundary, see remarks in [lUj. 
Moveover, for each T > 0, we denote by Qt, St, Qt the set Vt x (0, T),dVtx (0, T) 
and Vt X [0, T), respectively. The following Lemma gives the existence of strong trace 
of any entropy solutions to degenerate parabolic equations. 



NEUMANN BOUNDARY CONDITION FOR DEGENERATE EQUATIONS 



5 



Lemma 2.1. Let Assumption \2.2\ hold. Suppose that u G L°°{Qt) o-nd B(u) G 
L2([0,r];Fi(fi)) satisfying 

dtU + V -fiu) = AB{u) inV'iQT) 

and 

dtv{u) + V -qiu) - Ap{u) <0 inV'iQT), 
for any entropy- entropy flux triple (?7,q, where q' = 'r]'i',p' = rj'B'. Then u 
possesses a strong trace on the boundary. 

Remark 2.3. We notice that under the nonlinearity-difjusivity condition, the proof 
of this lemma is quite the same with the proof of Theorem 1.1 in [T7j. Indeed, under 
condition (12. 2p . the flux function f is genuinely nonlinear on the interval where the 
diffusion function B is degenerate. Thus, combining the technique used in [17j and 
utilizing the results of Vasseur [23], we can easily show that u possesses a strong 
trace on the boundary. 

2.2. Boundary-layer sequence Q and Normal Traces. In order to deal with 
the boundary terms, we introduce the concept of boundary-layer sequence, see [21] . 

Definition 2.3. We call {Cs} a boundary-layer sequence if Q ^ nC^{Q) is a 

sequence of functions such that 

lim C5(x) = 1 Vx G f]; < C<5 < 1; Cs = on dQ. 

5-^0+ 

Remark 2.4. For any boundary- layer sequence (s, —'^Cs converges to the outward 
normal n of the boundary, i.e., for any if) G {H^{VL))^ , 

lim I ^l)- VCs = hm / (V ■ ij)Cs = - I V = - I ifj ■ n. 

'5^0+ in ^^0+^ Jn J an 

Now we introduce the notion of divergence measure fields given by Chen and 
Frid in [iniin], see also [21]. Set 

VM\Qt) ■■= {F G L2(Qr,R^+i) : 3C> :| / F ■ {dt,V)<j)dxdt \< C||0|Uoo(2^)} 

JQt 

the space of L?'{Qt) vector fields whose divergence is a bounded Radon measure in 
Qt- Properties of divergence measure fields have been well studied by Chen and 
Frid in [inilll]. In particular, the following lemma which generalize the Gauss-Green 
formula for divergence measure fields holds, (see also [22]): 

Lemma 2.2. Let 7 G (L2(Qr))^+i be such that divj" is a bounded Radon measure 
on Qrp. Then there exists a linear functional 7^ on W2''^{Tjt) H C(S7^) which rep- 
resents the normal trace T ■ ^ on J^t in the sense that the following Gauss-Green 
formula holds: for all ip G C^{Qt), 

<%,i!>= [ ^(v-j^)+ / vv^-j^, 

Jqt J Qt 

and < 7^,-0 > depends only on -^Isy. 
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Let {Q} be a boundary-layer sequence, according to discussions in [22] we know 
that for any i/j € C^(Qt) 

< 7;, ^ >= - lim / :F- VCs^dxdt, (2.5) 



£-S>0 



T 



where G (L^(Qt))^ is the space part of J-" = (J-q, ^) and the hnear functional 7^ 
represents the normal trace JF ■ n. 

3. Uniqueness of Entropy Solutions 

In this section, we shall first give a definition of entropy solutions of problem 
( II. ip which is mainly involved with boundary condition (13. 5p . Then we show several 
useful lemmas which are mainly concerned with the techniques in dealing with the 
boundary terms appeared in the proof of Theorem 13. 1[ Finally we can prove its 
uniqueness. 

Definition 3.1. A function u G L°°(Qt) is called an entropy solution of the initial- 
boundary value problem (II. ip if the following conditions are satisfied : 

(1) (regularity) The following regularity properties hold: 

B{u) e L\0,T;H\n)) (3.1) 

VA;GM: {\u-k\,sgn{u-k){i{u)-f{k))-V\B{u)-B{k)\) eVM^Qr). (3.2) 

(2) (interior entropy condition) \/k G M, V0 G C^{Qt) with > 0, 

[ \u-k\dt(l)+{sgn{u-k){f{u)-f{k))-V\B{u)-B{k)\} ■V(pdxdt>0; (3.3) 

Jqt 

(3) (initial condition) The initial condition is satisfied as a limit in the following 



sense: 



ess lim / |it(x, t) — Mo(x)|dx = 0; (3.4) 
Jn 

(4) (boundary condition) For any boundary- layer sequence {^5} and nonnegative 
function if G C^(Qt), Vfc G R, 

lim / sgn(w - k){VB{u) - f(n)) ■ Vx^v^dxdt > 0. (3.5) 
Jqt 

Remark 3.1. By the theory of divergence-measure field, the entropy- entropy flux 
pair 

{\u - k\,sgn{u - k) (f (n) - f (A:) - V{B{u) - B{k)))) 
has normal traces for any A; G M, i.e., the limit 

lim / sgn{u- k){i{u) -i{k) -V{B{u) - B{k))) ■V^^s^d:>cdt (3.6) 

^-^0 Jqt 
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exists. On the other hand, by Lemma \2.1l the limit 

lim / sgn(M - k)i{k) ■ Vx^5<^dxdf (3.7) 

exists. Therefore, from (13. 6p and (13.71) . the limit in (13. 5 p exists. 
Remark 3.2. Taking k = ±||-u||loo in (13. Sp . we g'e^ 

lim / (f (m) - VB{u)) ■ Vx^V'dxdt = 0. (3.8) 

/t seems more natural to impose the condition (13. 8p as t/ie boundary condition. But 
it turns out to be too weak to ensure the uniqueness of entropy solutions (see [6jj. 
On the other hand, if the function {{u) — VB{u) has strong trace on the boundary, 
the equality (13. Sp reduces into the classical boundary condition 

{f (u) -VB{u)) -11 = on on. (3.9) 

Remark 3.3. If the problem (II. ip is fully degenerate, i.e, B'{u) = 0, the definition 
above is exactly the same as the definition of entropy solutions for zero-flux problem 
of conservation laws in j5]. 

Remark 3.4. It is easy to show that (13. 3p . (13. 4p and (13. 5p are equivalent to the 
following inequality: V/c G R,\fif G C^{Qt) with if>0. 



+ 



[ \u- k\dt^ + {&gn{u - k){i{u) - i{k)) - V\B{u) - 5(A:)|} Vy^dxdt 

Jqt 

I |mo(x) - fc|(^o(x)dx + I [ sgn(ii^ - A;)f (A;) ■ mpdn^-^dt > 0, (3.10) 
Jq Jo JdQ 

where v^o(x) = V^(0,x) and vT' is the strong trace of u on the boundary. In fact, 
(13.1 op implies (13. 3 p and (13. 4p obviously. Take (p = {1 — with i/j G C^(Qr) in 
(I3.10p and let 6^0, the inequality (13. 5p follows immediately. On the other hand, 
taking (f) = Cs(f in (13. 3p where ip G C^{Qt) o-nd using boundary condition (13. Sp . we 
get (I3.10p immediately. 

The following two lemmas are useful in the proof of uniqueness of entropy solu- 
tions. 

Lemma 3.1. Let u be an entropy solution defined in Definition l3A\ and D C 
be an open set such that D fl dfl = {x = {x,xn)\xn = 0} and D (1 Q = = 
{x, xn)\xn < 0}. Then for any boundary-layer sequence {Q} o,nd nonnegative func- 
tion i/j G C^{D X (0,T)), WA{x,t) G L°^{D H dQ x (0,T)), the following inequality 
holds: 

liminf / sgn(M - A{i,t)){y B{u) - i{u)) ■ VOV^dxdt > 0. (3.11) 
Jqt 
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Proof. First we claim that for any simple function S{x,t) = '^Eix{o,T)ki where 



ki G M, U™ = Dndn, and 

l£;iX(0,T) = 
the following inequality holds: 



i=l 



1, {x,t) e E,x {0,T) 
0, otherwise 



[ \u- S{x,t)\dtifi+ {sgn{u- S{x,t)){f{u) -i{S{x,t)) -VB{u))}\/ipdxdt 

JO.T 



T 

N-1, 



+ / sgn{u^ - S{x,t))i{S{x,t)) -mpdW^-'dty^S.U) 
Jo Jan 

where (f G C^{D x (0, T)). In fact, we only have to show that for each i, 
I I \u-ki\dtip+ {sg-n{u-ki){i{u) -i{ki) -VB{u))]Vipdyidt 

JxN J Ei 

+ I [ sgn{u^ -ki)i{ki) -ULpdn^-^dtyO. (3.13) 

Jo J Ei 

This follows easily by taking k = ki and ip = lEii^ in fl3.10p . Now, suppose A{x, t) G 
L°°(D ndflx (0, T)), then there exists a sequence of simple functions {An{x, t)}^^ 
such that 

\im An{x,t) = A{x,t) a.e. for {x,t) e D n dn x {0,T). 
For each An, we know that fl3.12p holds. Let n — )■ oo, we immediately get 

/ \u - A{i,t)\dt<p + {sgn{u - A{x,t)){i{u) - f{A{5:,t)) -VB{u))}V<pd:>cdt 
Jqt 

[ [ sgn(M^ - A{x, f))f t)) ■ n^dn^-Mt > 0(^.14) 
Jo Jdn 

Taking ip= {I- Cs)^^ in (1314]) where G x (0, T)) and letting 5 ^ 0, we get 

(13. lip immediately. □ 

The next lemma comes from and is modified in our case as follows: 

Lemma 3.2. Let D be a bounded subset in and 001,002, Ji, J2 ^ Let 
Pi{z) = ipiiz) where p G C^([— 1, 1], M) is a nonnegative even function such that 
J^p{r)dr = 1. Then 

limsup / / sgn(a;i(x) - uj2{y)){Ji{x) - J2{y))pi{x - y)dxdy 

i^oo JJdxD 



'Qt 

+ 



< / sgn(a;i - i02){Ji - J2)dx + / | Ji - J2|dx. 

' D J {a;i=tiJ2} 
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Moreover, if Ji = J2 a.e. on uji = uj2, then there exists 



lim // sgn.{ui{x)-U2{y)){Ji{x)-J2{y)pi{x-y)dxdy = / sgn(a;i-a;2)(^i-^2)cix. 
'^"^JJdxd Jd 

(3.15) 

Now we can show the uniqueness of the entropy solutions to problem (11.11) . 

Theorem 3.1. Suppose that u{'x,t) and i;(y, s) are entropy solutions of initial 
boundary value problem (II. ip with initial data uq and vq, respectively. Then for 
a.e. t > 0, 



u(x, t) — t'(x, t)|dx < / |iio(x) — t>o(x)|dx. (3.16) 
n Jn 

In particular there exists at most one entropy solution to the initial boundary value 
problem (II. ip . 

Proof. Let (^(x, t,y, s) be nonnegative smooth function which satisfies, 
e(-,-,y,s)eCo^(QT) for fixed (y,s)eQT, 

e(x, t, ; ■) G C^iQr) for fixed (x, t) G Qt- 
By standard "doubling of variables" (see ^), we get for test function ^(x, t,y,s) 

{\u{^,t)-v{y,sMt + Cs) 
+sgn(ti(x, t) - v{y, s)){f{u{^, t)) - f{v{y, s))) ■ (V.^ + VyO 

-(Vx|S(w)(x, t) - Biv){y, s)\ + Vy\B{u){^, t) - Biv){y, s)\) ■ (V^C + VyO}d^ > 0, 

where dz = dxdtdyds. 

For 2 = 1,2, let e C;f^(M^), < A < ^ with supp^i CC D, where D C is 
an open set such that 

D n (9^2 = {x = (x, xn)\xn = 0} 

and 

D n i7 = {x = {x,xn)\xn < 0}. 

We also take ip2 = ^ on the support of i/ji such that '0i(x)'02(x) = '?/'i(x). We 
denote ?/^i(x)?/^2(y) by V(x,y). Choose 9{t) G C^{0,T),e > 0. Let ^(x,t,y,s) = 
Csi^)Cri(j)pi{t - s)pm{x - y)pn{xN " ?/Ar)6'(t)^(x, y), where 0, Cr? are boundary-layer 
sequences and pi,Pm,Pn are sequences of mollifiers and x = {x,XN),y = {y,yN) are 
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local coordinates induced by i/j. Then we get 



Jj {\U- v\CsCvPl,m,n1pO' + Sgn{u - v){i{u) -i{v)) ■ V^+y1pC,5C,nPl,m,nd 
QtxQt 

-{V^\B{u) - B{v)\ + Vy\B{u) - B{v)\) ■ V^+yiljCsC,pi,m,ne}dz 



JJ Sgn(M - V){f{u) - i{v) - V^B{u)) ■ V ^C,5CnPl,Tn,nd'>P<^Z 
rxSr 

J J sgn(M - v){i{u) -i{v) + VyB{u)) ■ V yCr,Cs Pi,m,nOi^dz 



> 



QtxQt 

+ JJ V^\B{u)-B{v)\-Vy(r^CsPl,m,nOijdz 
QtxQt 



+ JJ Vy\B{u)~ B{v)\-V^CsCrjPl,m,nOi^dz 

QtxQt 
=: Ji + J2 + J3 + J4, 

where we denote piPmPn by pi,m,n for simplicity. Let 



Jl = Jl + J J Sgn{u - V)f{v) ■V^CsCr,OlpPl,m,ndz, 
QtxQt 

where 

Jl = J J sgn{u - v)(Vx-B(m) - f(n)) ■ V^CsCvdi'pi,m,ndz. 

QtxQt 

For a.e. (y, s) E Qt, let 

Fs{y, s) = sgn(M(x,t) - v{y, s)){V^B{u) - f(u)) ■ V^CsCr,^i^pi,m,nd^dt. 

J Qt 

By entropy boundary condition fl3.5l) . we have for a.e. (y, s) G Qt, 

limF5(y,s) >0. (3.17) 

(5->-0 

So, by Fatou's lemma, we get 

liminfJi> / liminf F5(y, s)dyds > 0. (3.18) 

Therefore, we conclude that 

liminf Ji > — / / sgn{u^ — v)i{v) ■ UxOippi mndV.^~^dtdyds, 
'?.'5^o J J 
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where Ex is the boundary domain with respect to x, nx is the outward normal vector 
to Sx and u'^ is strong trace of u on the boundary. Similarly, we obtain 

liminf J2 > / / sgn(w — v'^)i{u) ■ UyOi/jpi mndx.dtd'H^'^ds, 

^'V^O J J 

lim J3 = — / / sgn(M — v'^)\/xB{u) ■ n„6'?/^pi ^ ndxdM7/^~"^ds, 

^'V^O J J 

lim J4 = / / sgn(w^ - v)VyB{v) ■ n^9'iljpi.m,ndydsd'H^~'^dt. 

Thus we get 

jj {\U - v\pipmPntp0' + Sgn{u - v){f{u) - i{v)) ■ V^+y^pi,m,nO 
QtxQt 

-(Vxl^H - B{v)\ + Vy\B{u) - B{v)\) ■ V^+yi:pi,m,n9}dz 

> J5 + Jq, 

where 

J5 = - sgn(M - v^)iV^Biu) - f(M)) • nyetPpi^m^ndxdtdn^-^ds, 

and 

J6= jj sgn{u^ -v)iVyBiv)-iiv))-iiMPi,m,ndydsdn''-'dt. 
Next let /, m — )■ CX3, using lemma [3^ we obtain 



liminf J5 = — limsup / / sgn(tt — f'^)(Vx-B('u) — ^iu)) ■ iiyOiJj pi^^^nd^dtdT-L ds, 

l,m^oo l,m^oo J J 

> / sgn(ii - f^)(Vx5(M) - f(M)) ■ (0, -l)eiJiPnixN)dxdxNdt 

Jqt 

- [ \{V.,B{u) - f{u)) ■ (0, -l)9iJiPn{xN)\dS:dxNdt 

JE 

:= J7 + Js, 

where E := {{x,xi^)\u{x,xjyi) = v'^{x)}. Let a;„ = 2 J°^p„(^)d^, then we know 
VuJn = 2(0, — l)p„,(xAr). Therefore, using lemma ISTTl we get 

liminf J7 = - liminf / sgn(M — v'^)(V-B(w) — f('u)) ■ VuJnOipidxdxN > 0. 

n— ^-oo 2 n— >oo J Qj, 
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Next, we shall use the special property of rectangle domain to show that 

liminf Jg = 0. 

n—^oo 

We first note that 

V,B{u)-{O,-l) = -d,^B{u) = on E. 
On the other hand, we can easily show that 

liminf /" |/^|^p„(x^)^/;idx = 0. (3.19) 
Je 

Indeed, from remark [321 VV' E C^{{E\dn) n D), we have 

liminf / f n Pn{x n)'4' = ^ , 



n— >oo 



which implies /at = a.e. on dE fl dVL fl D. Thus, fl3.19p holds. Based on the above 
calculations, we get 

liminf liminf Js > 0. (3.20) 

Similarly, we have liminf liminf Jg ^ 0. Collecting all the limits above, we obtain 
the following inequality: 

/ \u- f |V'i(x)6'' + sgn(M - v){{{u) - {{v)) ■ V^iJiO 

JQt 

-(Vx|S(m) - B{v)\) ■ Vx^/'iMxdt > 0. (3.21) 

Now, we shall do the same partition of VL as in [16]. Let Vt = ^ '^t) be an 

open bounded rectangle with 2N faces 

{dn)i* = {{xi, ■ ■ ■ ,Xi_i,a*,Xi+i, ■ ■ ■ ,Xd);aJ < xj < a+ for j = 1,2, ■ ■ ■ , d, j ^ i} 

and the outward normal rij* to Q along {dQ)i* for i G {1,2,- ■■ ,d}, where the 
super-index * denotes the symbol + or — . We set Sj. = (0,T) x {dQ)i*. Set 
J = {1+, ■ ■ ■ , A^+, 1", ■ ■ ■ , A^-} and Jq = {0} U J. For > and i* e J we set U^,, 
{dfl)'^,, Q^*, Q^, and A^* as follows: U^* is the open subset of all x E fl such that 
dist{x,{dQ)i*) < u and dist{x, {dQ)i*) < dist{x, {dQ)j*) for all j G {1,2,-- - ,A^} 
with j 7^ i. {dQ)i* is the subset of all x G {dQ)i* such that x — srij. G U^, for all 
s G (0,z^). fi^. = {x — snj*;x G (9fi)^, ,s G (0,z^)}, the largest cylinder generated 
by n^. included in U.K. 9.^. = {x- sni.;x G (9fi)^,,s G A^. = Uji\n'^,. 

Thus, we have meas(Uj*gjAj',) < Const. z/^. Moreover, we set z* = if z = 0, 

Q^Q = Q\Ui*ejU^* and = Ui-gj^fiJ'.. Since the family {^7o^ fiJV, is an open 

cover of Q"^", we can choose {(/>i}{o<i<2Ar+i} a partition of the unity subordinated to 
the open cover. It's easy to see that for any j G {1, 2, ■ ■ ■ , 2N + 1}, supp^j fl dfl = 
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{x = {x,xn)\xn = 0} after change of coordinates. Using (13.211) with 0j in place of 
ijji and summing on i we get 

/ /_|M(x,t) - v(x,t)|^'dxdt > 0, (3.22) 

Jo JQ.'^'' 

for any 6 G C^(0, T). Let v goes to zero, we get 

/ |M(x,t) - t;(x,t)|0'dxdt > 0, (3.23) 

JQt 

Thus inequahty (13.161) follows from (I3.23P in a standard way. □ 

4. Existence of Entropy Solutions 

In this section, we prove the existence of entropy solutions for specially chosen f 
and B. Specifically, we assume that 

5'(0 = 0, for < ^ < Uc] B\^) > 0, for Uc<C<M. (4.1) 

and 

f(0 = 0, for e>Wc. (4.2) 

Remark 4.1. Note that conditions (14. ip and (14. 2 p are compatible with the nonlinearity- 
diffusivity condition (12. 2p . 

We shall use vanishing viscosity method to prove the existence result. For e > 0, 
we consider the following regularized parabolic problem: 

' dtu' + V ■ f (m^) = AB{u') + eAu', (x, t) e Qt, 
M^(x,0) = M§(x), ^eQ, (4.3) 
(f (m^) - VB{u') - eVu') ■ n = 0, (x, t) G S^, 

where is a sequence of smooth functions converging to Uq in L^{Q) assuming 
values in [e, M — e] for e > sufficiently small. The existence and uniqueness of 
classical solutions to (14.30 follow from standard arguments, see [[H], Chapter V]. 
Before we prove the existence of entropy solutions, we need a uniform a priori bound 
for u^. 

Lemma 4.1. Suppose that Assumption \2.1\ holds, then the approximate solutions 
satisfy < < M . Furthermore, there exists a subsequence {Sn} and a function 
u{-K,t) G L°°(Qt) such that u"" u in LI^^{Qt) as En and <u < M. 

Proof. Multiplying (14. 3p ^ by (-u^ — M)+ and integrating over Qt, using the boundary 
condition (14.30 0. we get 



(u'-M)+|Mx+ / V(5(ti) - S(Af)) ■ V(m - M)+dxdf 
n Jqt 

+e [ |V(M'-M)+pdxdt= /" f (m) • V(m' - M)+dxdt. 
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By Assumption 12.11 and using Young inequalities, we easily get 



{u' - M)+|Mx <C \{u'- M)+|Mxdt. 
n Jqt 

By Gronwall's inequality and noting that UttolU"" — ^? S^^ — 
same approach, using ("u^)" as test function, we get > 0. Therefore, by the 
compactness result of Lions, Perthame and Tadmor [20], we may extract a subse- 
quence of solutions of (14. 3p which converges in LI^^{Qt) to a function m(x, and 
0<u(x,t)<M. □ 

Lemma 4.2. The limit function u of solutions tt^ of the regularized problem fl4.3p 
satisfies ( 13. ip and (13. 2p m Definition \3. 1\ 

Proof. Multiplying (14. 3p ^ by B{u^) and integrating over Q^, we get 
I dtu'B{u')dii.dt+ ! (V ■ f (M"))S(M^)dxdt 

Jqt J Qt 

= / A5(m-') • S(u")dxdt + / 5Au^S(ti^)dxdt. (4.4) 

Jqt J Qt 

Integrating by part and using (l4.3P o. we get 

J Qt JQt 

£M2, 



|V5(w')rdxdt+ / eWu' ■WB{u')dii.dt. (4.5) 

Qt JQt 

Let A be primitive function of B, i.e.. A' = B. Then by Lemma 14. 1\ we get 

- [ |VS(ii^)pdxdt + e [ B'{u')\Vu'\'^dxdt 

< [ dtA{u')dxdt + - [ |f(M=)|2dxdt < C, 
Jqt 2 Jqj, 

where C is a constant independent of e. Thus (13. 2 p holds. To show the stated 
VAi^ property, we use the same approach as in [21j (see also |6]). First define 
approximation sgn^ and | ■ |r; of the sign and modulus functions by 

sgn^(^) := \ ^> kl < V, ■= / sgn^(C)dC, ^ > 0. 

1-1, rK-r], ° 
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Multiply dSDi by sgn^{B%u^) - B%k))ip where k e R and ip E C^(Qt) and 
integrating by parts over Qt, using ( ]4.3P n. we get 

Qt 



Qt 



{i{u') - f{k) - V{B%u') - B%k))}sgn^{B%u') - B%k))Vijdxdt 



+ I {i{u')-i{k))sgn'^{B'{u')-B'{k))V{B%u')-B'{k))i)dxdt (4.6) 

JQt 

-I \u' - k\r,tp \^ dx+ [ {u' - k)sgn'^iB%u') - B' {k))dt{B' {u') - B'{k))ijdxdt 

Jq J Qt 

+ f {u' - k)sgn^{B'{u') - B'{k))dtipdxdt 

JQt 

+ [ [ f{k) ■ nsgn^{B%u') - B'{k))^dn'' 
Jq Jan 



'0 Jan 
ji + /2 _^ ^ j6 



Now we consider the limit of the righthand side of (14.61) for 77 — )■ 0. First we note 
that sgn('u^ — k) = sgn(i?^(M^) — B^{k)) due to the monotonicity of -B^(-) and by 
using Lebesgue's theorem, we get as — )■ 0, 

II [ {sgn{u' - k){f{u') - f(A;)) - V\B%tf) - B%k)\}Vijdxdt. 

JQt 

To estimate we shall use the fact that usgn'^(M) < X{u:0<\u\<r]} (see [21]) and 
noting that sup|^|<;y^ B'^'iO — ^5 "^^ get 

- i{k)>gn'^{B%un - B%k))V{B%u) - B%k))\ 
< '""^^^^-f^^^ ViB^iu^) - B%k))\xxie,,) 
where X[e,rf) = {{x,t) : < \B'^{u'^) — B'^{k)\ < 77}. Therefore, 

|^2| < sup|g|<Mf(0 ^|^^|^^^^^^ r |v(s^(„^) _ B%k))\dxdt. 

£ Jl{e,-n) 

Observe that the meas X(e, ?]) — )■ as — )■ 0, so — )■ as 77 — )■ 0. Next, we see 
that as — )■ 0, 

/3^/3:= / \u'{x,T)-k\'ilj{x,T)-\u'Q{x,T)-k\^{x,0)dx 
Jn 
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and thus |/q| < 2(M+ ^) ||'C||l°°(Qt)- The integrand of satisfying 
= \iu' - k)sgn'^{u' - k)dt{B%u') - B%kM 

< IdtiB'iu") - S^(A;))X{(x,i):0<|«-(x,t)-fe|<r,}||V'l|L-(QT)- 

The same argument employed for imphes that — )■ as — t- 0. Finally, using 
Lebesgue's theorem, we get 



J5 ^ J5 := /" \u' - k\dtil)^xAt 
Jqt 



and 



rp 

I^^I^-= I f f{k) ■nsgniB'iu') - B'{u))^d'H^. 
Jo Jdn 

It's easy to see that |Jo| < T|(9fi||f(A;)|||'0||2,oo(gy). Collecting the estimates on to 
yields that all terms of the righthand side of (14 .Op possess a limit as — )■ and 
are in particular uniformly bounded with respect to 77. Therefore, taking -0 = 1, we 
see that there exists a constant Ci, depending possibly on e (but not on 77), such 
that 

[ {VB%u') - B%k) fsgn'^{B%u') - B%k))dxdt < Ciie). 

JQt 

Consequently, the sequence 



{i?,,J,>o := {{VB^{u^) - B^{k)YsgT^^{B%u^) - i?^(fc))},>o 

is bounded in V-{Qt) with respect to t] and therefore also in C(Qt)', the dual space 
of C{Qt) of continuous functions on Qt- By compactness of the weak-:*r topology of 
C(Qt)', we deduce that, up to subsequence, the sequence {i?e,r?} converges towards 
an element G C{Qt)' in the weak-^ topology. Thus for any ip G C°°(Qt) we can 
pass to the limit 77 —t- in (14. 6 p to get 

<Ee,^p> = I {sgn{u' -k){i{u^) - f{k)) -V\B'{u') - B'{k)\}Vijdxdt 

JQt 

- I {\u'{x, T) - k\tlj{x, T) - \u'{x, 0) - fc|0(x, 0)}dx (4.7) 



e 



k\dtilJ+ [ [ i{k)-nsgn^{B'{u')-B'{k))ipdH^. 
Jo Jdn 



+ / \u' 

'Qt Jo Jdn 
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On the other hand, since sgn^ > 0, we have E^^n > for every 5,77 > 0. Therefore, 
we get 

\<E„'ilj>\ 
M\l^{Qt) 

= hm — ^ [ (VB^u') - B%k)fsgn'JB%u') - B%k))4)dxdt 

v^o \\iP\\l--(Qt) Jqt 

< hmsup f {VB%u') - B%k))hgn'^{B%u') - B%k))dxdt 
Thus, we get from (14. 7p with ip = 1 

<Ee,ij_>_ ^ f ll^.(^^y) _ \u%x,0)-k\}dx 



+ [ [ f{k) ■ nsgn^{B%u') - B%k))ijdn^ (4.8) 
Jo Jan 

Using estimate from Lemma (I4.ip . we deduce that there exists a constant C2 which 
does not depend on e such that | < i?e,'(/' > I < C'2||'0||Loo(gy) for all e > 0. 
Consequently, E,, is bounded in C^Qt)' and up to subsequence E^ converges in 
the weak-^ topology to a function E G C(Qt)'- Now we pass to the limit e — )■ 
in (14. 7p . Note that — k\ converges strongly to \u — k\ in C{0,T] L^{Q)) and 
sgn{u'^ — k){f{u'^) — f (A;) — V-B'^(ii^)) converges weakly in L'^{Qt) to sgn{u — k){{{u) — 
i{k) - VbIu)). Thus we conclude that for all (p e C^{Qt), 

<E,(j)>= I {\u-k\dT(f> + sgn{u-k){i{u) -i{k) -VB{u)) ■V(t)}dxdt (4.9) 

JQt 

Since i? is a Radon measure, we obtain from (14. 9 p that for all G C^{Qt)-i 



{\u-k\dt<p + sgn{u-k){i{u) -i{k) - W B{u)) - W <i)}dxdt < C||0||loc(q^). (4.10) 

'St 

This in particular implies the VAi^ property (13. 2p . □ 

The existence of entropy solutions is presented in the following theorem. 

Theorem 4.1. Suppose 'Uo(x) ^ [0, Af] for a.e. x G and Assumption \2. 1\ and \2. S\ 
hold. Moreover, we require that assumptions (14. ip and (14. 2 p are satisfied. Then the 
limit function u{x,t) in Lemma 4-i is an entropy solution of initial-boundary value 
problem (II. ip . 

Proof We divide our proof into two steps. 

Stepl. Let ri G C^iR) be a convex function and 0(x, t) G C^([0, T)xQ). Multiplying 
(gD, by77'(«^)0(x,t),weget 



77(M^)0tdxdt+ / (q(w^)-Vp(M^) + eVr7(M^))- V0dxdt+ / 77(wg)0(O, x)dx > 0, 
Jqt Jn 

(4.11) 
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where q' = 77'f' and jo' = rj'B'. Let 5 — )■ 0, we get 



?7(u)0tdxdt+ / (q(ti) - Vp(u)) ■ V0dxdt + / r/(Mo)0(O, x)dx > 0. (4.12) 
St ■''St 

By approximation, the above inequahty still holds for r]{u) = \u — k\ for any k eM.. 
So, conditions (2) and (3) hold. 

Step 2. Let (/?(x,t) e C^(Qt)- Multiplying (jODi by v?(x, t) and sgn(M^ - A;)(^(x, t) 
respectively and integrating over Qt, using fl4.3P n and (14. ip . we get 

u'dtipdxdt+ [ {i{u')-f{k)-VB{u')-eVu'}-Vipdxdt = 0, (4.13) 

St -^St 
and for k > Uc 

I \u'-k\dtipdxdt+ [ sgn{u'-k){i{u')-i{k)-VB{u')-eVu'}-Vipdxdt>0. 

J St "'St 

(4.14) 

Let e — )■ 0, we have 

/ udt^dyidt+ {i{u) -f{k) -VB{u)}V(fd^dt = 0, (4.15) 

JSt "'St 

and for k > 

u - k\dtipd:>cdt + / sgn(w - A:){f (m) - f (A;) - VS(n)} ■ Vy^dxdt > 0. (4.16) 

St >^St 

Take f{x,t) = (1 - C<5(x))V^(x, t) in (ITOj) and (Hl6|) with ?/^(x,t) G Co~(Qt) and 
then let 5 — )■ 0, we have 

lim /" (f(w) - V5(w)) ■ VO^dxdt = 0, (4.17) 
is.T 



and for k > Ur 



lim / sgn(M - k){VB{u) - i{u)) ■ V^V^dxdt > 0. (4.18) 

•^-^o isT 

Next we shall show that (KTE\) holds for k < u^. Using (H^]), (H^]) and fHTTjl . we 
have for k < Uc 

lim/ sgn(n - A;)(f(M) - V5(ii)) ■ VO^/' 



i5-s>0 



{u>Uc}U{u<Uc} 



lim / (f (m) - VB{u)) ■ VCsip + lim / sgn(M - k)i{u) ■ VCsifj 



= lim / (sgn(^z -k)- l)i{u) ■ VCs^ = 0. 

The last equality in the above formula holds since VB{u) = on the set {u < Uc} 
and thus f('U^) ■ n = a.e. on the set d{u < Uc} fl dQ by (14.170 . Combining 
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the regularity results from Lemma 14.21 the existence of entropy solutions follows 
immediately. □ 
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